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The slow compaction of a gently tapped granular packing is reminiscent of the low-temperature 
dynamics of structural and spin glasses. Here, I probe the dynamical spectrum of granular com- 
paction by measuring a complex (frequency-dependent) volumetric susceptibility \ v . While the 
packing density p displays glass-like slow relaxations (aging) and history-dependence (memory) at 
low tapping amplitudes, the susceptibility \ v displays very weak aging effects, and its spectrum 
shows no sign of a rapidly growing timescale. These features place x v in sharp contrast to its di- 
electric and magnetic counterparts in structural and spin glasses; instead, \ v bears close similarities 
to the complex specific heat of spin glasses. This, I suggest, indicates the glass-like dynamics in 
granular compaction are governed by statistically rare relaxation processes that become increasingly 
separated in timescale from the typical relaxations of the system. Finally, I examine the effect of 
finite system size on the spectrum of compaction dynamics. Starting from the ansatz that low fre- 
quency processes correspond to large scale particle rearrangements, I suggest the observed finite size 
effects are consistent with the suppression of large-scale collective rearrangements in small systems. 

PACS numbers: 45.70.Cc, 64.70.P- 



RANDOMLY PACKED GRAINS AS A MODEL 
GLASS 

Filling a container by pouring in spherical grains nearly 
always leads to a disordered packing. The "ground state" 
of the system, a face-centered cubic packing and the dens- 
est arrangement of spheres, is rarely achieved unless put 
in by hand. Shaking or tapping the container, causing 
the contents to settle and compact, works only so well: 
the packing remains disordered, and its packing density 
p approaches but does not exceed the random close pack- 
ing (RCP) density p rcp sa 0.64 [flUEI]. Despite being a 
robust state that is readily reproducible, little is known 
about the origins of random closed packing. 

Surprisingly, the detailed arrangement of particles in 
RCP is remarkably similar to the molecular arrangements 
in a simple liquid H ED E0 BD E] ; but unlike a liquid, RCP 
is mechanically rigid and resists shear — it is jammed 
[lOj . The same combination of rigidity and disorder also 
characterizes structural glasses, such as window glass and 
hard candies. Like RCP, a glass is an elastic solid whose 
molecular arrangement is nearly indistinguishable from 
that of a liquid. Indeed, the structure of simple metallic 
glasses are well-modelled by RCP [HI US US]. Draw- 
ing on these similarities, Liu and Nagel proposed that 
jammed packing of particles and glass- formation in sim- 
ple liquids may be placed into a unified Jamming Phase 
Diagram spanned by inverse density 1/p, temperature 
T, and shear stress (Fig. [T} [14. The region of high 
density, low temperature, and small shear is occupied by 
jammed or glassy states; the rest are filled by fluid states. 
Jamming is the consolidation of loose grains into a rigid 
packing when p is increased at T = 0, whereas the glass 
transition is the freezing of a liquid into a glassy solid at 
a finite T . The principal aim of this work is to explore 
the relationship between jammed grains and glasses by 
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FIG. 1: (Color online) The Jamming Phase Diagram spanned 
by orthogonal axes inverse density 1/p, temperature T, and 
shear stress o xy . The shaded volume is occupied by jammed 
or glassy states; the rest is filled by fluid states. Point J marks 
the T = 0, p — p rC p jamming transition. 



adopting methods and principles from glass experiments 
to study granular materials, in an effort to identify shared 
features, as well as differences and new effects. 

It is attractive to think random closed packing and 
glasses, unsolved problems on their own, may in fact be 
closely connected. It also invites skepticism to think sys- 
tems so different microscopically can exhibit close analo- 
gies. We are therefore encouraged to find that granular 
packing is not only glass-like in structure, it can also ex- 
hibit slow relaxations reminiscent of the low-T dynam- 
ics of structural glasses, and of other "glassy" systems 
such as spin glasses and kinctically constrained computer 
models [15j [TBJ [17]. For glass- forming liquids such as 
glycerol or viscous Si02, small reductions in T can lead 
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FIG. 2: (Color online) Granular compaction under steady 
tapping; v = 1/p is the specific packing volume, t is the 
number of taps applied, and T is the tapping amplitude (see 
Chapter 2.1). Strongly tapped (F large), the packing quickly 
reaches a steady state; weakly tapped (r small), the packing 
volume compacts slowly, and the steady state is not accessible 
within the experimental timescale. 



to large (and apparently diverging) increases in the re- 
laxation time (and the fluid viscosity). Below the glass 
transition temperature T g , the characteristic relaxation 
time of the liquid becomes greater than the experimental 
timescale: the liquid falls out of equilibrium and becomes 
a glass. Because it is out of equilibrium, the properties of 
the glassy state are not stationary but depend on the ex- 
perimental timescale (aging) and on the thermal history 
of the system (memory). For a packing of macroscopic 
particles with strongly dissipative interactions (friction), 
thermal energy is unimportant, and the dynamical evo- 
lution of the system is driven by external energy input 
(such as tapping or shaking). Strongly shaken, a granu- 
lar packing quickly reaches a loosely-packed steady state. 
Weakly tapped, the packing compacts slowly in a log- 
arithmic fashion with time, and the steady state may 
not be accessible on an experimental timescale (Fig. |2j; 
furthermore, the detailed evolution depends on how the 
packing was prepared PH EU US [HJ |22J [23] . 

One may therefore take the random packing of grains 
- despite its athermal, dissipative nature — to be a 
model "glass" , in terms both of structure and of dynam- 
ics. This raises the question whether there are shared 
physical principles connecting the various microscopically 
disparate systems commonly regarded as "glassy". Sim- 
ply put: is a granular packing glassy in the same sense 
structural and/or spin glasses are glassy? If there is a 
deep connection, then granular packings may be a con- 
venient experimental system for investigating the origins 
of glassy behavior. Here, I have developed a volumetric 
spectroscopy, modeled after the dielectric and magnetic 



spectroscopies of structural and spin glasses, to directly 
probe the spectrum of relaxation processes in granular 
compaction, with goal of ultimately comparing it, in de- 
tail, with the spectrum of glasses. Because granular pack- 
ing is an athermal system without a well-defined statis- 
tical mechanics, analogies with the behavior of glasses 
at finite T cannot be justified a priori. I will adopt a 
purely operational approach: analogies will be proposed 
naively, and they will be justified solely on the light they 
may shed on the experimental data. 

This paper is organized as follows: 

Section II describes the volumetric spectroscopy and 
gives an operational definition for the "susceptibility" \ v . 
I will then describe how the spectroscopy is performed, 
and Xv measured, in a granular compaction experiment. 

Section III describes the behavior of \ v as a function of 
the taping amplitude (which is the principal control pa- 
rameter for granular compaction) and the frequency (the 
dynamical spectrum). Despite some similarities to the 
dielectric/magnetic susceptibilities of structural and spin 
glasses, Xv is distinctive in exhibiting very weak aging. 
More surprisingly, there is no indication in the spectrum 
of Xv of a diverging relaxation time as the tapping am- 
plitude is turned down. Comparing the compaction pro- 
cess with the dynamical spectrum, one finds that granu- 
lar compaction is characterized by a broad separation of 
timescales. 

Section IV suggest that the apparently distinctive be- 
havior of Xv may be generic to a complex configurational 
specific heat. In analogy to the well-known relation- 
ship between specific heat and energy fluctuations, the 
relationship between Xv and specific volume fluctuation 
is found to be consistent with an effective Fluctuation- 
Dissipation Theorem. Furthermore, computer simula- 
tions of two standard spin glass models found the com- 
plex specific heat in these models appear to mirror the 
behavior of Xv- If Xv can indeed be taken as a "spe- 
cific heat" , it suggest the slow dynamics of granular com- 
paction is governed by statistically rare relaxation pro- 
cesses which become increasing decoupled from the typ- 
ical dynamics of the system. 

Section V describes how the system size affects the 
dynamics of granular compaction. For small systems, 
Xv is suppressed within an intermediate frequency band 
and within an intermediate window of tapping amplitude; 
this is consistent with the suppression of large-scale cor- 
related motions in small systems. But for sufficiently 
weak driving, the low-frequency responses of small sys- 
tems are enhanced; I suggest this is the consequence of 
particle/wall interactions at the boundary. 

To conclude, in Section VI, I will summarize the most 
salient results and briefly discuss their implications for 
understanding "glassy" dynamics. I will close by sug- 
gesting ways that this work can be extended to address 
several important open questions. 
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THE VOLUMETRIC SPECTROSCOPY OF 
GRANULAR COMPACTION 

The dynamics of structural and/or spin glasses is com- 
monly studied via an appropriate spectroscopy. Here, 
temperature T is usually the principal parameter that 
controls whether the system is in the glassy state. Hav- 
ing set the temperature, the dynamics of system is char- 
acterized by some susceptibility, which describes how the 
system responds to small perturbations. For example, 
the dielectric susceptibility e, often used to characterize 
glass- forming liquids, describes the response of the net 
polarization P when a small electric field E is applied: 
e = dP/dE. If the perturbation can be applied in a 
frequency-controlled way, e.g. with an AC electric field, 
a dynamical spectrum of response versus frequency can 
be obtained. Such a spectrum will encode the relative 
distribution of relaxation times in the system. 

In a similar spirit, I will characterize granular com- 
paction using a complex (frequency-dependent) volumet- 
ric "susceptibility" \ v , which describes the response of 
the specific packing volume v = 1/p to small changes in 
the external driving (tapping), and extract the dynami- 
cal spectrum of granular compaction. This should allow 
us to make a detailed comparison between the slow com- 
paction of a granular packing and the low-temperature 
dynamics of structural and spin glasses. 

An operational definition of \ v 

In a typical compaction experiment, a granular packing 
(held inside a container) is compacted with a sequence of 
identical taps, separated sufficiently apart so that the 
packing comes to a complete rest before being tapped 
again. Each tap is one single period of sinusoidal (30 Hz, 
in my case) vertical vibration, and the tapping amplitude 
r is given by the peak-to-peak acceleration a pp of the 
tap, normalized by g — 9.8 m/s 2 : T = a pp /g. In order to 
extract a susceptibility, the packing must be perturbed at 
the same time it is being tapped; however, experimentally 
there is no readily accessible secondary field which can 
couple to the system. 

The solution here is to use tapping both as the prin- 
cipal control parameter and as the perturbing "field" . 
Instead of a sequence of taps with identical amplitudes, 
the tapping amplitude is modulated about a mean level 
[Fig Ha)], 

r(t) = r + (5rsin(27r/i). (1) 

The mean amplitude Tq serves as the principal control 
parameter, while the modulation amplitude ST serves as 
the perturbing "field"; / is the modulation frequency in 
inverse taps, and t is time measured in the number of taps 
applied. In response to the modulated tapping, the tem- 
poral evolution of the specific packing volume v = l/p 



can be described as the sum of a slowly varying com- 
ponent vq, and a frequency- locked component which is 
phase-shifted with respect to the input [Fig |3jb)], 

v(t) rnv + 5v sin(2ir ft- cp). (2) 

Here 5v is the response magnitude, and <p is the phase 
shift between the modulated tapping and the volumetric 
response. One can then define the real and imaginary 
(in- and out-of-phase) parts of a complex susceptibility 

Xv {XviXv) as 

Xv = (jfj cos Lp, (3) 

x'l = (^) sin ^ W 

This susceptibility \ v describes how the specific volume v 
responds to a small change in the tapping amplitude <5r. 
For sufficiently weak perturbations, the response should 
be linear, 5v oc ST, and Xv should be independent of the 
modulation aamplitude. 

In practice, the phase shift ip and the response am- 
plitude Sv are not separately extracted from the v(t) 
data. Using the orthonormality of trigonometric func- 
tions, Xv can be directly extracted by projecting v(t) onto 
its in- and out-of-phase components (with respect to the 
tapping modulation). Once a v(t) time series has been 
collected, it is divided into time-blocks of length 1//, 
synchronized to the tapping modulation. Within each 
time-block, v(t) is discretely transformed to extract Xv 
directly: 

xU*) = 3f E v(t + T)sm(2irfr), (5) 

T=0 

x"W = S E w(* + r)8ia(27r/r-7r/2). (6) 

T=0 

Formally, Eqs.[5][6]are exact projections only if the back- 
ground vq of v(t) is stationary. However, since vq is 
almost always slowly- varying on the timescale of 1//, 
Eqs. [5] [6] are excellent approximations. 

Experimental design 

The experimental design here is partly modeled after 
that of rcf [18 . I use a cylindrical packing whose height is 
much larger than its diameter in an effort suppress large 
scale convections during the compaction process. The 
cell is a smooth-walled cylindrical tube made of poly- 
carbonate, 50 cm tall, with an inner diameter D = 2.5 
cm. For experiments, the cell is filled with particles to a 
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FIG. 3: (Color online) (a) Amplitude modulation of discrete taps, (b) Input: T(t), amplitude-modulated tapping as a function 
of the tap number t (solid line). Experimentally measured response: frequency-locked modulation in the specific packing volume 
v(t) (filled circles); here the time series has been divided into two time-blocks with lengths of 1/f — 50 taps. Note that F(t) 
and v(t) are slightly out of phase. 
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FIG. 4: (a) The experimental setup used to measure x v 
(schematic not drawn to scale), (b) A photograph of the pack- 
ing; its image profile is fitted to extract the packing height, 
and hence the specific packing volume v. 

height of ps 40 cm [Fig. |4^a)] . To minimize static build- 
up, the inside cell wall is coated with a hard conduc- 
tive polymer coating (Agfa Orgacon S103, ITO America 
Corp., Tempe, AZ). The cell is mounted on a computer- 
controlled electromechanical shaker, which is carefully 
leveled and set on a vibration-damping platform. To 
minimize unwanted lateral motions during tapping, the 
experimental cell is braced to a stiff external framework 
with heavy-duty rubber bands. 

The particles used are spherical zirconium oxide/silica 
ceramic (material density = 3.88 g/cm ) particles (Glenn 
Mills Inc., Clifton, NJ). Compared to more commonly 
used glass beads, these are mechanically much tougher. 
They also have a much smoother surface, which mini- 
mizes abrasion on the cell and on the particles them- 



selves. These features help to ensure good reproducibil- 
ity over the course of very long experiments and millions 
of taps. Particles of three different diameters were used: 
d = 0.5 mm, 1.1 mm, and 1.8 mm, all with ~ 10% poly- 
dispersity. Because air effects can significantly alter the 
compaction process 24J, particularly when particles are 
small, once loaded the cell is sealed and evacuated to a 
pressure < 5 Pa. 

To measure the specific packing volume, a photograph 
of the packing is taken after each tap using a computer- 
controlled digital camera running Astro IIDC (Outcast 
Software, Calgary, Alberta). The cell is backlit with a 
uniform, bright light source. Where the cell is filled with 
particles (which are opaque), the image is completely 
dark. By fitting the intensity profile of the image to 
an error function, the location of the top surface of the 
packing can be found [Fig. |4jb)]. This is then compared 
with known reference marks on the cell to extract the 
packing height, and hence, the packing volume. The ab- 
solute accuracy in the determination of v is limited by the 
unevenness is the top surface and by the accuracy with 
which reference marks can be placed. These are compen- 
sated by the large fill height, which keeps the absolute 
uncertainty in v to ~ 5 x 10~ 3 . The relative accuracy, 
which is what matters most here (since the signal I seek 
consists of small changes about a slowly-varying back- 
ground), is somewhat better (~ 10 -4 ) due to the fact 
that image profiles can be tracked with sub-pixel reso- 
lution. Compared to methods that measure the packing 
density locally, such as capacitance [18] or 7-ray absorp- 
tion |20j . the imaging technique measures the bulk av- 
erage and therefore displays much less fluctuation from 
tap to tap. This is necessary in order to obtain a clean 
measurement of Xv ■ 

To ensure a fresh start at for each new experiment, the 
packing is "reset" by tapping at T = 21 for at least 10 3 
taps. The result is a reproducible loose packing that has 
achieved a steady-state, where all memory of its prior 
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tapping history appears to be lost. Unless noted other- 
wise, this is the loose-packed initial state I will use for 
each new experimental run. 



CHARACTERIZING THE VOLUMETRIC 
RESPONSE 

The volumetric susceptibility \ v depends on both Tq, 
which is the principal control parameter, and on the mod- 
ulation frequency /. It should not depend on ST if one 
is in the linear response regime, i.e. Sv oc ST. At low Tq, 
where granular compaction exhibits glassy behavior, Xv 
should in principle depend on the experimental timescale 
(aging), and the previous history of the packing (mem- 
ory) . Indeed this is what we expect based on the behav- 
ior of dielectric and magnetic susceptibilities in structural 
and spin glasses. Unless explicitly noted, all the results 
described in this Section were obtained using d = 0.5 mm 
particles, giving an effective system size of D/d = 50. 



Dependence of \ v on To and ST 

To measure Xv as a function of T (at a fixed frequency 
/), the mean tapping amplitude is ramped from To = 21 
to Tq = 2 and cycled back in steps of Aro = 0.5 and a 
dwell time per step of Td w = 500 taps. Fig.[5]plots x„(r ) 
measured at different frequencies, where the modulation 
amplitude is set to ST = 1 . Because little hysteresis is ap- 
parent during the cycling, I have averaged together the 
values of \v obtained during the down-ramp with that 
obtained during the up-ramp. As a function of Tq, both 
x' v and x'v vary non-monotonically with Tq, exhibiting a 
broad maximum at an intermediate range of Tq for all 
the frequencies measured. This shows intermediate Tq is 
where compaction is most sensitive to variations in the 
tapping amplitude. The broad maximum is structured, 
with sub-peaks at To = 10 and Tq = 7.5, and a "knee" 
at To = 5. These sub-features are found in both x'v and 
x", and are located at the same values of T . At lower 
frequencies, one obtains a stronger the real response, and 
the fine features in Xvfto) stand out more pronouncedly. 
On the other hand, the imaginary response at low fre- 
quencies is quite weak. 

As the mean tapping amplitude Tq and the tempera- 
ture T are, respectively, the principal control parameters 
for granular compaction and for glass formation, it is 
reasonable to compare X«(r ) with the T-dependent sus- 
ceptibilities of glasses. Qualitatively, Xv(^o) has broad 
similarities to the magnetic susceptibility xm(T) of spin 
glasses, such as CuMn, as well as the dielectric sus- 
ceptibility e(T) of structural glasses such as glycerol 
[231 |2"51 |2"5] . However, for most structural and spin 
glasses, the peak in e(T) or xm{T) has a simple struc- 
ture, whereas Xv(Xo) exhibits multiple sub-peaks. For 
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FIG. 5: (Color online) The (a) real and (b) imaginary parts of 
the volumetric susceptibility \v as a function of To, measured 
at different modulation frequencies /. 



glasses, the peak in e(T) or xm(T) shifts to lower T for 
lower frequencies. In granular compaction, the fine fea- 
tures X«(ro) remains at the same values of Tq. 

Fig. [6] compares Xv(^o) measured using ST = 1 with 
that obtained using a smaller modulation amplitude ST = 
0.5; the results are practically identical. This indicates 
the experiments are largely in the linear response regime 
within the experimental range of / and Tq. On closer 
examination, the sub-features in Xv(Po) are more pro- 
nounced at the smaller ST. This suggests the sub-peaks 
of Xv mark values of To where the packing is particularly 
sensitive, and where the response is most nonlinear. This 
is reminiscent of the extra sensitivity one expects of sys- 
tems near a critical point. A similar effect can be found 
in spin glasses, where low-magnetic field (weaker per- 
turbation) measurements yield a taller, sharper peak in 
the magnetic susceptibility. But where in the spin glass, 
the peak of xm shifts to higher T when measured using 
smaller magnetic fields |30l 131] , in granular compaction 
the sub-peaks in Xv do not shift with ST. 

For structural or spin glasses, where the dielectric 
or magnetic susceptibility peaks is often taken to mark 
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FIG. 6: (Color online) Testing for linear response: (a) x'v and 
(b) x" measured using ST = 1 (filled symbols) and SF — 0.5 
(open symbols) at several different frequencies (as indicated 
in figure legend). The similarity in the resulting curves shows 
the experiment is in the linear response regime. 



the glass transition temperature T g . Here, we can ask 
whether the features in x v correspond to any "transi- 
tions" in compaction dynamics with Tq. 

A candidate location for such a transition is the branch 
point r*, where the reversible and irreversible branches 
of the volumetric curve v(Tq) converge [19] . As shown 
in Fig. [7] starting from a loose-packed initial state and 
ramping up Tq slowly from a small value, the pack- 
ing initially compacts with increasing Tq, before it re- 
verse course and begins to dilate with further increases 
of Tq; this trajectory is the "irreversible branch" . If Tq is 
then slowly ramped down from large amplitude to small, 
the packing compacts monotonically, becoming denser at 
lower Tq; subsequent cycling of Tq between strong and 
weak tapping (at the same ramping rate) will remain 
on this so-called reversible branch. For Tq > T* , the two 
branches converge, and the packing reaches a steady state 
with a unique volume on the experimental timescale (~ 
inverse ramping rate). For To < T*, w depends not only 
on Tq but also on prior tapping history. In this sense, T* 
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FIG. 7: The volumetric curve v(To), showing the irreversible 
(open symbols) and reversible (filled symbols) branches. Ar- 
rows indicate the direction To is ramped, and T* marks where 
the two branches appears to converge. 



is like T g in that for T < T g , the properties of a glass de- 
pends on its thermal history. As the convergence of the 
irreversible and reversible branches is gradual, T* is a 
not precise location but depends (weakly) on the rate at 
which r is ramped [T5] . This is reminiscent of the fuzzi- 
ness of T g , which depends weakly on the rate at which 
the glass-forming system is cooled down from high T. 

The susceptibility Xv(Tq) can be thought as the com- 
plex "derivative" of the reversible branch of the com- 
paction curve v(Tq). Thus the broadly peaked form of 
Xvfto) is not a surprise given the form of w(Fo). The 
sub-peak in \v at Tq = 10 is located near, but slightly 
below r*. But for the equally prominent peak located at 
Tq = 7.5, there is no outstanding feature in w(Fo) that 
appears to correspond to it. 



Aging and memory effects in x„ 

Fig. [8ja) shows the behavior of granular compaction 
subjected to a step- wise shift protocol in Tq, starting 
from a loosely-packed initial state. With the initial onset 
of tapping at r = 7, the specific volume v exhibits a 
strong aging signature and slowly compacts in a quasi- 
logarithmic fashion with t. After 10 4 taps, Tq is turned 
down to Tq = 5, and the system responds with renewed 
bout of aging and relaxes to a smaller volume. After 10 4 
taps at r = 5, the system is returned to r = 7, where- 
upon the packing dilates to a volume close to what it 
had just before the excursion to Tq = 5. When the two 
To = 7 portions of the data are stitched together, the 
result is identical (aside from a transient) to an uninter- 
rupted r = 7 aging curve. 

In structural and spin glasses, the experiment de- 
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FIG. 9: (Color online) The relaxation of x'v (solid lines), x'v 
(broken lines) at T = 7 for / = 10" 1 (orange), / = 100" 1 
(purple) , and / = 1000" 1 taps" 1 (black). Xv is mostly re- 
laxed after a timescale on the order of 1//. 



timescale ft, starting from loosely packed initial states; 
here T n = 7, and / = 10"\l0" 2 , and 10" 3 taps" 1 . 
The susceptibility appears to be largely relaxed after a 
timescale on the order of 1//. 



FIG. 8: (Color online) Granular compaction subjected to an 
step-wise shift protocol ro = 7 — ► 5 — ► 7, starting from a 
loose-packed initial state, (a) Specific volume v(t) (crosses), 
showing strong aging signatures both initially, and at the neg- 
ative shift to To = 5. When the two To = 7 portions of v(t) 
stitched together (grey-filled circles), the results is the same 
as an uninterrupted Fo = 7 v(t) reference curve (bold broken 
line), aside from a brief transient, (b) The volumetric sus- 
ceptibility x v measured in the same experiment. The aging 
signature is very weak, and the memory effect is trivial. 



scribed above is carried out via step-wise shifts in the 
temperature T, and the aging and memory effects are 
observed in the dielectric or spin susceptibility . 
It is a classic demonstration of both the aging and mem- 
ory effects in a glassy system. For granular compaction, 
in the same experiment where aging and memory effects 
are so clearly displayed in the specific volume v, the vol- 
umetric susceptibility x v exhibits no discernible aging 
signature [Fig.[8^b)]. Xv relaxes within a few modulation 
periods 1 // to apparently static values in response to 
the shifts in IV In this respect, Xv behaves rather unlike 
dielectric and magnetic susceptibilities of structural and 
spin glasses with their strong aging effects. A closer ex- 
amination shows that aging effects in \v are not entirely 
absent, but that they are very weak. Fig. [9] plots the re- 
laxation of Xv at different frequencies against the reduced 



The spectrum of x„ 

The weak aging effects in \v can be employed to our 
advantage in the measurement of its spectrum. Even if 
the system is not in steady state, the measured suscepti- 
bility, due to its weak aging effects, should be close to the 
"true" steady-state value. Thus the spectrum of Xv can 
be measured down to fairly small To, where the steady 
state is experimentally inaccessible, with the confidence 
that the resulting spectrum will be close to the steady- 
state one. Here, starting from a loosely packed initial 
state, T is slowly ramped stepwise (step size Ar = 0.5, 
dwell time/step Td w = 500 taps) to the working Tq. The 
packing is then tapped at the working Tq (without modu- 
lation) for an extended period until cither steady state is 
achieved, or until the volumetric evolution is slow on the 
order of the lowest frequency to be probed. Only then is 
the amplitude modulation (ST — 1) turned on, and the 
frequency swept from high frequencies to low. 



Fig. 10 plots the dynamical spectrum Xv(f) obtained 
at different values of Tq. Here, the imaginary part of 
the spectrum x"if) measures the spectral distribution 
of relaxation processes in granular compaction. In over- 
all form, the spectrum of Xv is similar to its dielectric 
and magnetic counterparts in glasses, including a peaked 
imaginary part. As with glasses, the spectral peak in 
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FIG. 10: (Color online) The (a) real and (b) imaginary parts 
of the Xv spectrum, obtained at different values of To (as in- 
dicated in figure legend). 



Xv is wider than that of a single Debye peak. In the 
vicinity of the peak frequency x'v ^ can be fitted to the 
generalized Havriliak-Negami form [33] 



Xv(f) 



"X.V 



A.-v 



(7) 



[1 + (^Tt/t)"]/ 3 ' 

with a 1=3 1 and « 0.5 (note a — /3 — 1 recovers 
the Debye form) [Fig. 11 'a)]; however, the low frequency 
wing of x'v is i n excess of what a generalized, single-peak 
form (such as Eq. [7]) can account. For very small To, the 
response of v(t) to low- frequency amplitude- modulated 
tapping can be strongly anharmonic. This can result in 
anomalous features in the Xv spectrum for small To at 
low /, such as negative values of x'v- 

It is clear that while the spectral peak in x"(f) varies 
in height and shape with To, the peak location is nearly 
stationary. This is in sharp contrast to the dielectric and 
magnetic spectra of glasses, where a small reduction in T 
can shift the spectral features to dramatically lower fre- 
quencies [551 157]. By fitting the spectral peak in x" 
to Eq. |TJ the peak frequency / p (and equivalently, a char- 
acteristic relaxation time of the system r p = l// p ) can 



FIG. 11: (Color online) (a) The spectral peak in x" a t To = 7. 
Lines are fits to the spectral peak using Debye (broken line) 
and Havriliak-Negami (solid line) forms, (b) The character- 
istic relaxation time t p = l// p plotted against 1 /To- 



be extracted. Fig. [TT|b) plots the inverse peak frequency 
l// p = T p versus l/r obtained over a broad range of 
r . For r > 7.5, t p w 25 taps does not change with r ; 
below To = 7.5, t p grows modestly from t p w 25 to 75 
taps, and appears to saturate for very low Tq. There is 
no indication that t p diverges as Tq — * 0. 

The characteristic relaxation time of t p w 25 — 75 
taps can only be described as remarkably short. Even 
at the largest Tq within this range, an loosely-packed ini- 
tial state will take thousands of taps to reach a steady 
state. Below To w T* w 11.5, when the steady state 
is practically inaccessible, the characteristic timescalc t p 
extracted from x'v(f) is at least 3 orders of magnitude 
smaller than the time it takes for a loosely-packed initial 
state reach the steady state. The behavior of t p with To 
in granular compaction is in sharp contrast to the analo- 
gous behavior of glasses, inverse peak frequency of the di- 
electric and magnetic susceptibilities diverge as t p ~ e x l T 
or even faster. 
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An anomalous "susceptibility" 



The behavior of Xv, insofar as it describes the responses 
of a system with "glassy" dynamics, is highly anomalous. 
The weakness of its aging behavior, as well as the ab- 
sence of a diverging timescale in its spectrum, render \v 
sharply distinct from the familiar dielectric and magnetic 
susceptibilities of structural and spin glasses. Strangely, 
Xv does not appear to capture the glassy dynamics of 
granular compaction: there is a broad separation of scale 
between the characteristic timescale, extracted from the 
dynamical spectrum, and the time it takes loosely-packed 
initial state to arrive at steady state. 

It may be that the anomalous behavior of Xv is a 
unique feature of granular compaction, and has no coun- 
terparts in structural or spin glasses glasses. In ther- 
mal systems, a susceptibility is strictly defined by the re- 
sponse of a macroscopic thermodynamic observable with 
respect to its conjugate control parameter or field, e.g. 
Xm = dM/dh. For granular compaction, where there 
is no well-defined statistical mechanics, no principle dic- 
tates Xv should be the analog of any thermodynamic sus- 
ceptibility. But as I will suggest in Sec. Ill, the distinc- 
tive features of Xv may not be so anomalous after all; 
they may instead represent the generic behavior of the 
configuration specific heat in glassy systems. 



VOLUMETRIC SUSCEPTIBILITY AS A 
"SPECIFIC HEAT" 



Granular compaction is an athermal system with 
strongly dissipative dynamics; nevertheless, I suggest 
tapping plays the role in granular compaction akin to 
temperature in glasses. Most importantly, tapping acti- 
vates transitions from one static configuration of particles 
to another, thereby allowing the system to explore con- 
figuration space; this is reminiscent of the temperature- 
driven exploration of the phase space in thermal systems. 
Additionally, in granular compaction, the tapping ampli- 
tude r is the principal control parameter that deter- 
mines whether the dynamics is glassy, while in glasses 
the analogous role is taken by T. 

With this in mind, I propose the following analogy: 
taking To as a temperature-like control parameter, with 
v playing the role of energy density as proposed by Ed- 
wards et al. [35j [36] , then Xv should thought as the 
complex "specific heat" of granular compaction. While 
this analogy cannot be justified from first principles, op- 
erationally it proves to be very helpful towards under- 
standing the experimental results. 
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FIG. 12: (Color online) Testing the effective FDT relationship 
between volume fluctuations and the volumetric susceptibility 
at To = 11.5. Plotted are the power spectrum of steady-state 
volume fluctuations S(f) (light gray crosses), S(f) smoothed 
by averaging over logarithmically spaced bins (bold broken 
line), and the combination x"(f)/f (filled symbols). 



Connection to specific volume fluctuations 

In thermal equilibrium, the specific heat describes the 
energy (or equivalently, entropy) fluctuations of the sys- 
tem via the Fluctuation-Dissipation Theorem (FDT). 
Even though a tapped granular packing is not an equi- 
librium thermal system, it is reasonable to suppose On- 
sager's Regression Hypothesis — that the system does 
not distinguish between spontaneous fluctuations and 
those induced by a small external perturbation — should 
still hold to some approximation, particularly as the ex- 
periments described here are largely in the linear re- 
sponse regime [57]. If Xv is a specific heat-like quantity, 
with v playing the role of energy density, then Xv should 
also describe steady-state volume fluctuation in granular 
compaction. Modeling after the FDT between specific 
heat and energy fluctuations [35J [5^3 > I hypothesize the 
steady state of granular compaction can be described by 
an effective FDT 



S(f) = NT^x"U)lf, 



(8) 



where S(f) is the power spectrum of specific volume fluc- 
tuations, TV is the number of particles in the system, and 
T e g — T e g (/, Tq) is the effective temperature that charac- 
terizes the steady-state fluctuations on a timescale 1//. 
In the simplest case, in exact analogy to FDT in ther- 
mal equilibrium, T c ff should be frequency-independent 
for / < 1/tq, where r is a microscopic relaxation time. 

The experimental set-up here is not well-suited for 
measuring volume fluctuations. The bulk- averaging over 
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a large packing of small particles, necessary in order to 
yield a clean measurement of Xv , also averages out fluc- 
tuations in v(t). In addition, the steady state is only 
accessible for large r , where x" is small and subject to 
large uncertainties. The range of Tq over which all these 
conflicting demands can be met, and over which Eq.[8jcan 
be tested, is accordingly very narrow. With these caveats 
in mind, I have independently measured steady-state vol- 
ume fluctuations and the susceptibility spectrum Xv(f) 
at r = 11.5. Fig. 12 plots the power spectrum S(f) 



of specific volume fluctuations; plotted on top of S(f) is 
the combination x"(/)//- While S(f) and x"U)lf are 
clearly distinct at very high frequencies, they are roughly 
consistent with each other at lower frequencies over 2 
decades in bandwidth. This suggests away from the very 
high frequencies, the relationship between x.v and spe- 
cific volume fluctuations is consistent with that of an ef- 
fective FDT, with a T e g that is frequency-independent. 
The high-frequency discrepancy between S(f) and Xv/f 
could be an artifact caused by tap-to-tap variations in 
the topography of the top surface, which register as spu- 
rious fluctuations in the pack height (and hence specific 
volume). It may represent an unavoidable departure as 
/ — > 1 taps -1 , due to the discrete nature of time in gran- 
ular compaction. 

I note here the S(f) shown in Fig. 12 has a different 
slope from that obtained by the earlier study of Nowak et 
al. [TS]. Here S(f) ~ f~ l , whereas ref. [2T] measuring at 
a similar tapping amplitude, found S(f) ~ /~ 2 over the 
same frequency range 0D]. The reason for this difference 
is not clear, but it could due to finite size effects. The 
system measured in Ref. [19] consists of d = 2.2 mm glass 
particles confined in a D — 1.8 cm cylinder (D/d — 8.2), 
where the measurement here is of d — 0.5 mm ceramic 
particles confined in a D — 2.5 cm cylinder (D/d = 50). 



Comparison with glass-forming liquids 

While the specific is not usually thought as a dynami- 
cal quantity, the complex specific heat (at constant pres- 
sure) c p — (cf p , c'p) has been measured experimentally for 
a number of glass-forming organic liquids j-JU 0U 321 HI] ■ 
Approaching T g from above, the real part c' p (T) is nearly 
flat until close to T g , when it drops to a lower value over a 
relatively small span in temperature. Accompanying the 
stepwise drop in c' p (T), the imaginary part c^(T) grows 
from zero far above T g to a peak centered at T g , be- 
fore decaying to zero with further cooling. The spectral 
evolution of c p with T shows a peak frequency in c^(/) 
that shifts rapidly to lower frequencies (indicative of ever 
larger relaxation times) as T — > T g from above. Quali- 
tatively, the complex specific heat c p (f) of glass-forming 
liquids is similar to that of dielectric or magnetic suscep- 
tibilities; it does not behave like the volumetric suscepti- 
bility Xv(f) °f granular compaction. 



It is possible there exists a low-frequency band in the 
spectrum of x„, outside of the experimental window, 
which does slow down as Tq — > 0. This low-/ band could 
then be analogous to the spectral features observed in 
the complex specific heat of glass- forming liquids, and 
would be responsible for the glassy character of granu- 
lar compaction. In a molecular dynamics study of liquid 
SiC>2 approaching T g from above, Scheidler et al. [45] 
found the spectrum of the constant-volume specific heat 
c"(/) consists of two bands, a high-/ band which does not 
shift with temperature, and a low-/ band which shifts 
to lower frequencies as T decreases. There, the high- 
/ band comes from the vibrational motion of molecules 
about a metastable configuration, while the low-/ band 
is contributed by configurational rearrangements. It may 
be that the spectral features observed in granular com- 
paction is analogous to the high-/ band of vibrational 
processes in liquids. But as granular compaction takes 
place in an athermal system with zero kinetic energy, it 
is unclear this analogy is appropriate. 



Comparison with spin glass models 

There is a crucial distinction between glass-forming liq- 
uids and granular compaction in the nature of their en- 
tropy. For a glass-forming liquid, its entropy has both 
configurational contributions (from the spatial arrange- 
ment of the molecules) and dynamic contributions (from 
the kinetic motion of the molecules). The entropy of 
granular compaction, on the other hand, is purely config- 
urational. Because specific heat is intimately connected 
to energy and entropy, this difference between glassy- 
forming liquids and granular compaction may prove cru- 
cial. Therefore it may be more appropriate then to com- 
pare Xv with the specific heat of purely configurational 
systems, such as spin glasses. 

The DC magnetic specific heat c mag (where the elec- 
tronic and phonon contributions to the specific heat have 
been removed) has been measured for several spin glass 
systems via scanning calorimetry. These have found 
c mag (T) to take broadly peaked forms, qualitatively sim- 
ilar to Xv(^o), which in some cases exhibits sub-peaks 
and "knees" [131171 HSJ. However, being DC measure- 
ments, these contain no dynamical information against 
which one can make a detailed comparison. 

In the absence of experimental measurements of the 
complex specific heat in real spin glasses, I have resorted 
to Monte Carlo computer simulations to study the purely 
configurational complex specific heat c = (c',c") in two 
standard spin glass models. I examined 2D Ising models 
governed by the Hamiltonian 

H = Y, JijSiSj -h^Si=H - KM, (9) 
(i,3) i 
where Si = ±1 are Ising spins occupying the sites of a 
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FIG. 13: (Color online) The imaginary part (c") of the spin 
glass specific heat £ spectrum, obtained at different temper- 
atures T for (a) the Edwards- Anderson (EA) model, (b) the 
±J model. The system size is 16 2 . Time is measured in Monte 
Carlo steps (MCS); 1 MCS = 256 spin flip attempts. 



16 x 16 square lattice with periodic boundary conditions. 
Here h is an externally applied magnetic field, and Jij 
are quenched random couplings between nearest neigh- 
bor spin pairs For the Edwards- Anderson (EA) 
model, {Jij} are distributed in two unit- variance gaus- 
sians centered at J = ±1; for the ± J model, takes on 
values of { — 1, 1} with equal probability. 

Fig. [13] plots the spectrum of the imaginary specific 
heat c"(/) at different temperatures for both the EA and 
±J model. For both models, the height and the shape 
of the spectral peak varies with T, but the inverse peak 
frequency l// p = r p grows only modestly and does not 
appear to diverge as T — > 0. Qualitatively, the overall be- 
havior of c"(f) in both models is similar to that of Xv(f)> 
and the forms of r p (l/T) in the spin glass models are re- 
markably close to T p (l/ro) as observed in granular com- 
paction (Fig. 14 1. I have also found at low temperatures, 
c shows a weaker aging effect than both the overall energy 
H and the magnetic susceptibility xm, reminiscent of the 
weak aging signature in the volumetric susceptibility Xv- 



FIG. 14: (Color online) The characteristic energy relaxation 
time t p = l// p , where / p is the peak frequency in c"(/), 
for the EA and ±J models (circles). Also shown are mag- 
netic relaxation times obtained from the magnetic suscepti- 
bility spectrum for the same models (diamonds). As T is 
decreased, there is a rapidly growing separation of timescale 
between energy (obtained from specific heat) and magnetic 
(from the magnetic susceptibility) relaxation times. 



Thus apparently distinctive features of Xv — weak ag- 
ing effects and the absence of a diverging timescale - 
also appears to characterize the complex specific heat of 
purely configurational spin glasses. 



Separation of timescales 

Like the volumetric susceptibility Xv in granular com- 
paction, the complex specific heat is characterized by 
short characteristic timescales, and does not appear to 
capture the glassy dynamics of the spin glass models 
studied. On the other hand, the magnetic susceptibil- 
ity in the same spin glass models are characterized by 
rapidly growing timescales that tend to diverge as T — > 



(Fig. 14 1. This shows on cooling, there is a rapid separa- 



tion of timescales between characteristic energy/entropy 
fluctuations, which remain relatively fast, and magnetic 
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fluctuations, which slow down enormously. However, all 
dynamical processes in the spin glass, including magnetic 
fluctuations, must ultimately contribute towards the spe- 
cific heat; that the spectrum of c shows no indication of 
a rapidly growing characteristic timescale suggests the 
glassy aspects in the low-T dynamics of spin glass must 
be governed by statistically rare processes. 

Is a similar effect at work in granular compaction? As 
Xv appears to describe volume fluctuations via an effec- 
tive FDT, then the fact that its spectrum shows no sign of 
a diverging timescale suggests the tapping-induced volu- 
metric compaction of granular packing cannot be under- 
stood on the basis of typical relaxations of typical volume 
fluctuations. If x,v can be thought as a specific heat, it 
should capture all the dynamics of granular compaction. 
That one sees no indication of a diverging characteris- 
tic timescale in x" suggests the glassy aspects of gran- 
ular compaction are controlled by statistically rare pro- 
cesses which become increasingly separated in timescale 
from typical relaxation processes as Tq is reduced, lead- 
ing eventually to the breakdown of crgodicity [HO [50] . 
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FIG. 15: The compaction curve wfTo), showing both irre- 
versible (open symbols) and reversible branches (filled sym- 
bols), obtained for packings with D/d = 50, 22.7, and 13.9. 
The ramping of To is in steps of Aro = 0.5, Td w = 500 taps. 



Rescaling the system: the effect on v and \ v 



THE EFFECT OF SYSTEM SIZE 



Specific heat is an intensive quantity that should be 
independent of system size. But when a system is con- 
fined with one or more very small dimensions, its be- 
havior can become qualitatively different from the bulk, 
and normally intensive quantities can become sensitive 
to the system size. Most obviously, the properties of a 
small system with large surface-to-volume ratio will be 
dominated by that of the interface, which may be qual- 
itatively different from that of the bulk. More subtly, if 
collective effects are important to the dynamics of the 
system, then a small system may be unable to support 
those collective modes characterized by large dynamical 
correlation lengths, and the contribution of such modes 
to the specific heat will be lost. 

Growing evidence from both simulations and experi- 
ments suggest that the collective motions of many parti- 
cles play a crucial role in the low temperature dynamics 
of structural glasses. In that case, finite size and confine- 
ment effects should modify the glassy state in nontrivial 
ways. As a model "glass" where the small system size is 
readily accessible, granular compaction is almost ideally 
suited for looking into how finite size effects can glassy 
behavior. The volumetric susceptibility Xv, as a specific- 
heat-like quantity, is then a potentially sensitive indicator 
of whether collective particle rearrangements are impor- 
tant to the compaction dynamics. 



By packing particles of different diameters d into a 
cylindrical cell whose diameter is fixed at D = 2.5 cm, 
the system is effectively rescaled as D — > D/d. Here, 
three different particles sizes were used: d = 0.5, 1.1, 
and 1.8 mm; this gives us packings with rescaled diame- 
ters D/d = 50, 22.7, and 13.9. Since the height (w 40 cm 
in all cases) of the packing is much larger than the cell di- 
ameter, D/d is the smallest system dimension and should 
be principally responsible for any finite size effects. 



Fig. 15 plots the volumetric curve i>(ro) for all three 



system sizes D/d, obtained using the same ramping, 
Ar /rdw = 0.5/500 taps. Aside from an overall shift to 
larger v (lower p), the f(Fo) curve for small D/d is not 
sharply different from that for large D jd. The reversible 
branch is slightly less steep, and the branch point T* is 
less clearly defined, but on the whole it is not obvious 
that any drastic changes have taken place. 

The differences stand out far more clearly when x v (Tq) 
for different system sizes D /d are compared against each 
other. As shown in Fig. [16] when the tapping is either 
very strong or very weak, x v (^o) appears to be inde- 
pendent of system size. But in the intermediate regime, 
where x v (Tq) rises to form a broad peak, finite-size effects 
manifest themselves as a systematic suppression of the 
overall feature, as well as the sub-peaks, with decreasing 
D/d. Over this intermediate range of To, a large system 
is more sensitive, and responds more strongly to small 
variation in the tapping amplitude, than a small system. 

Just as finite-size effects appear most prominently 
within an intermediate range of Tq, there is also a win- 
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FIG. 16: The effect of finite size on the Fo-dependence of \ v measured at frequencies of (a, b) / = 10 1 taps 1 , (c, d) / = 50 1 
taps -1 , and (e, f) / = 250 -1 taps -1 . The systems sizes are D/d = 50 (black), D/d — 22.7 (gray), and D/d — 13.9 (light gray). 



dow in frequency where it most strongly modifies the 
spectrum of Xv Fig. [17| compares the frequency spectra 
Xv(f) at fixed T for different system sizes D/d. For large 
r , the spectrum of Xv is essentially the same for all D/d 
[Fig. 17 a)]; but as r is lowered, the spectral peak in 
x'l becomes increasingly suppressed for smaller systems. 
The suppression is asymmetric: the low-frequency wing 
of the spectral peak is more strongly suppressed than 
the high frequency wing, so that for smaller D/d, the 
peak center is located at higher frequencies [Figs. 17 b, 



c)]. At still smaller values of To, while the spectral peak 
of a small system remains strongly suppressed, its low- 
frequency wing decays slower than that of the large sys- 
tem, so that at very low frequencies, x" f°r small D/d 
is systematically enhanced relative to that of large D/d 
[Figs.[l7^d, e)]. Finally, for very weak tapping, the spec- 
trum of Xv once again is essentially the same for all sys- 



tem sizes [Fig. 17 )] 



While the finite size effect manifests itself in the spe- 
cific volume as mostly a uniform shift to larger v (lower 
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density), its effects on \ v appears most prominently 
within windows both in Tq and in /. The smaller system 
has a faster characteristic relaxation time, as denoted by 
the inverse peak frequency of x"(/)> but given the non- 
trivial changes in the dynamical spectrum with D/d, it is 
difficult to state categorically whether small systems are 
more or less "glassy" . What is clear is that the packing 
with small D / d is not simply a shrunken version of the 
larger packing — a qualitative change in the dynamics 
has taken place. 



Extrinsic and intrinsic length scales 

For a granular packing confined inside a cylindrical cell, 
the surface-to-volume ratio scales with the extrinsically 
imposed cell diameter as (D/d) -1 . For the smaller sys- 
tem, those processes that take place near the boundary 
will exert a proportionally stronger influence on its over- 
all behavior. The wall not only restricts the range of mo- 
tion for particles near the boundary, it also restricts the 
ways these particles can pack together. This is a particu- 
larly strong effect when D/d is small and the curvature of 
the boundary proportionally large. This geometric hin- 
derance to efficient packing is likely the reason behind 
the reduced packing density when D/d is small. 

At low tapping amplitudes To < 7, it is observed that 
the particles next to the wall arrange themselves into an 
ordered triangular layer whose configuration appears to 
be nearly static from tap to tap. The larger surface-to- 
volume ratio of small systems means a greater proportion 
of the particles are semi-immobilized by the wall-induced 
ordering. The slow rearrangement of particles adjacent to 
the wall may explain the excess response of small systems 



at very low frequencies and weak tapping [Figs. 17 'd, e)] 



but this does not explain the asymmetric suppression of 
the spectral peak in x"(f)- 

Suppose the tapping-induced particle rearrangements 
during compaction are collective in nature and can be 
characterized by a set of correlation lengths {£} (which 
will depend on Tq). When the system size is small, those 
responses characterized by £ > D/d will be suppressed. 
Since it is not possible to track individual particle motion 
in these experiments, the presence of spatially correlated 
particle rearrangements cannot be directly ascertained; 
one must seek instead some indirect probe of the puta- 
tive length scales intrinsic to the system. Here I will rely 
on the following ansatz: at a given r , it is the collec- 
tive rearrangement of the particles on a large length scale 
that contributes to the low-/ responses, while small scale 
rearrangements give rise to the high-/ responses. In that 
case the correlation lengths can be mapped onto frequen- 
cies, £ — > £(/), and the spectral distribution of relaxation 
process, as measured by x"(/)> a l so describes the distri- 
bution of dynamic correlation lengths £. 

The asymmetric suppression of the spectral peak in 



x", when D/d is small, can then be understood as the 
suppression of large-scale collective rearrangements. As- 
suming that for a small system with strong finite size 
effects, the low-/ roll-off in the spectral peak takes place 
at / = fo, corresponding to when £(/d) = D/d, then a 
rough bound may be placed on the most probable cor- 
relation length £ p (corresponding to the peak frequency 
/ p ) for collective rearrangements in the largest system. 
From this admittedly crude bound, I estimate that for 
D/d = 50, £ p is not more than ~ 20 particle diameters. 
Since finite size effects are only apparent in an intermedi- 
ate range in To, £ must have a non-monotonic variation 
with To; £ must be smaller than the smallest D/d (=13.9) 
when the tapping is either very strong and very gentle, 
while reaching some maximum size at intermediate val- 
ues of Tq that mark the sub-peaks in Xu(r). 



Comparison with glasses 

While some aspects of the modified volumetric re- 
sponse in small systems can be attributed to the larger 
surface-to-volume ratio, the most prominent effect — the 
asymmetric suppression of the spectral peak in x'vi SU S~ 
gests a small system lacks some of the responses found in 
systems with a larger D/d, particularly if \ v is taken to 
be the "specific heat" of granular compaction. Presum- 
ably, these missing responses come from collective par- 
ticle rearrangements characterized by correlation lengths 
£ > D/d. From the variation in the finite size effect with 
tapping amplitude, these correlation lengths should have 
a non-monotonic To-dependence. £ should be small when 
the tapping is either very strong or very weak, and should 
grow to maximum size at some intermediate To ~ 10 
where the volumetric response is strongest. 

These conclusions agree qualitatively with what we 
know about dynamic correlations in supercooled liquids. 
Molecular dynamics simulations of supercooled liquids 
have found string-like clusters of mobile particles that re- 
arrange collectively; the characteristic sizes of these clus- 
ters grows from ~ 5 particle diameters far above T g to 
~ 20 particle diameters just above T g [HJ [52] . Experi- 
mental studies of glass-forming liquids using confinement 
within a nano-porous substrate [53J 153] . specific heat 
spectroscopy [55], multi-point susceptibilities [55], and 
NMR |57j all point indirectly towards growing dynamic 
correlations on the scale of tens of molecular diameters 
as T — > T g . Collective rearrangements have been directly 
observed in colloidal glasses where the glass transition is 
controlled by the colloidal volume fraction <fr. Approach- 
ing the glass transition point <j> g , the correlated clusters 
grow in size from a few particle diameters to a few tens of 
particle diameters, but the decrease rapidly as the system 
enters deeply into the glassy state (<p > 4> g )\EEl HH]- I n 
granular systems, string- like clusters of dynamically cor- 
related particles, on scale ~ 10 particle diameters, have 
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been directly observed in granular beds fluidized by gas 
flow; there the clusters grow in size as the external driv- 
ing is turned down 60J . Self-organized clusters have also 
been observed in 2D systems composed of particles that 
slowly swell in size and jam, but are otherwise not driven 
by external forces; here the clusters were found to grow 
in size as the system approaches the jammed state |61j . 



DISCUSSION AND OUTLOOK 



In this paper, I have described the experimental study 
of tapping-induced granular compaction, using a volu- 
metric spectroscopy to directly probe the distribution of 
relaxation processes. This extracts a complex "suscepti- 
bility" Xv, which describes how the specific packing vol- 
ume v = 1/p responds to small variations in the tapping 
amplitude T about a mean level. The principal goal of 
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this work is to make a detailed comparison between the 
dynamics of granular compaction with the slow relax- 
ations of structural and spin glasses. 

Glassy state in compaction dynamics 

Despite some qualitative similarities, the volumetric 
susceptibility \v is sharply distinct from its dielectric 
and magnetic counterparts in structural and spin glasses 
in important ways. \ v exhibits very weak signs of ag- 
ing, and its frequency spectrum shows no indication of a 
rapidly slowing relaxation time. There is a large separa- 
tion between te characteristic relaxation time extracted 
from Xv i an d the time it takes a loose-packed initial state 
to reach steady state. The slow, glassy dynamics of gran- 
ular compaction simply does not appear to manifest itself 
in the spectral response. 

Monte Carlo simulations of two standard spin glass 
models show that the complex specific heat c in these 
models bear striking similarities to x v . In particular, 
as T is reduced, the spectrum c shows only a modestly 
growing relaxation time in a manner that is remarkably 
similar to the behavior of Xv when T is turned down. This 
suggests the features of Xv are not distinctive to granular 
compaction, but may be generic to the specific heat in 
purely configurational glassy systems. This idea if further 
supported by the relationship between Xv and volume 
fluctuations, which appears to obey an effective FDT in 
analogy to the relationship between specific heat and en- 
ergy fluctuations in thermal systems. Interestingly, gran- 
ular compaction appears to split the difference between 
structural and spin glasses: structurally it is similar to 
former, while dynamically it behaves like the latter. 

If Xv can be taken as a specific heat, then it should con- 
tain contributions from all relaxation processes, including 
those responsible for glassy behavior. That a diverging 
timescale cannot be discerned from the spectrum of x v 
(or from the specific heat in spin glass models) suggests 
glassy behavior is largely controlled by a small minority 
of relaxation processes. As To or T is reduced, these rare 
glassy processes become rapidly separated in timescale 
from typical relaxations of the system, leading eventually 
to broken ergodicity and the departure from equilibrium. 

If my interpretation here of finite size effects in Xv 
are correct, it suggests that the glassy dynamics of com- 
paction is governed by rare events and characterized by 
a distribution of dynamical correlation lengths. These 
appear to point toward a "droplet" picture of the glassy 
state in granular compaction. For structural glasses, the 
droplet model suggests the glassy state is a mosaic struc- 
ture composed of small domains (droplets) drawn from a 
family of disordered equilibrium thermodynamic phases; 
the surface energy of the droplets scales weakly with 
droplet size, so that the droplets grow slowly and no sin- 
gle droplet can dominate the system 155] . Most of 



the relaxation processes in the system describe equilib- 
rium droplet fluctuations, and have a correlation length 
on the scale of droplet size. However, droplets also grow 
or shrink via rare, activated events, and these are the pro- 
cesses responsible for glassy dynamics. The results here 
suggest a similar droplet picture may be constructed for 
granular compaction as well. 

Connection to jamming 

A principal motivation of the present work is to explore 
the connection between structural glasses and granular 
packing, as suggested by the Jamming Phase Diagram. 
Jamming is the consolidation of loose particles into a 
rigid packing when the density is increased at T = 0. 
For frictionless hard spheres, there is a well-defined jam- 
ming point J (where p — p rcp ) that exhibits special scal- 
ing properties that mark it as an unconventional critical 
point [64, 65, 66 . Being a zero temperature phenomenon, 
the jamming transition does not inherently contain any 
dynamics. But evidence from experiments and simula- 
tions suggest that the critical properties of point J ex- 
tend beyond its immediate vicinity to influence dynamics 
at finite-T, where point J is not directly accessible, and 
be felt in a modified form as the glass transition [57]. 

In granular compaction, the ideal transition between 
jammed and unjammed states at point J is never ex- 
plored, since any packing one produces experimentally is 
automatically jammed. Because real particles have fric- 
tion, instead of a single jamming point at RCP, there is a 
large configuration space of jammed states, with densities 
ranging from random loose packing (RLP) at p r j p ss 0.55 
up to /9 rcp w 0.64. Tapping provides the dynamics for the 
system to explore this configuration space, with point J 
being an idealized limiting point. 

If the influence of point J can extend to finite-T and 
manifest itself as the glass transition, then a similar in- 
fluence may also be felt at finite friction and control 
the transition between glassy and steady state dynamics 
in granular compaction. I suggest the fuzzy transition 
between liquid and glass, and the fuzzy transition be- 
tween steady-state and glassy dynamics in granular com- 
paction, may be thought of as two aspects of the same 
idealized jamming transition, manifested non-ideally at 
finite T and finite friction, respectively. This may be 
why similar dynamical manifestations, such as collective 
rearrangements on the scales of few tens of particles, 
are found in both glasses and externally-driven granu- 
lar systems. In jamming, there appears to be a diverging 
dynamical correlation length associated with vibrational 
modes of the system [65] ES ES] . This diverging length 
scale may not survive way from zero temperature and 
zero friction, but its vestiges could be the source of the 
small-to-modest dynamic correlation lengths in glasses 
and in granular compaction. 
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Outstanding questions 

What are the rare processes and slow fluctuations that 
give us the glassy behavior of granular compaction? In 
spin glasses, it is the magnetic fluctuations that slow 
down on cooling. In structural glasses, viscosity diverges 
on cooling, leading to increasingly slow relaxation times 
for internal stresses. A good candidate for an analog in 
granular packing are fluctuations in contact force network 
between particle [70] |7T]. For very weak tapping, one ex- 
pects a contact network, once formed, will be minimally 
perturbed from tap to tap, and that the fluctuations in 
the contact forces will be characterized by very long re- 
laxation times. In this case, a susceptibility may not be 
available, since the obvious coupling field, gravity, is not 
normally controllable; but it should be possible to mea- 
sure the power spectrum of contact force fluctuations. 

What arc is role of kinetic degrees of freedom in the 
glass transition of liquids? Why does the complex spe- 
cific heat c p of glass-forming liquids show indications of 
a diverging timescale, while specific heat of spin glass 
simulations and the volumetric susceptibility x„of gran- 
ular compaction do not? One way this question may 
be attacked is to extend volumetric spectroscopy to con- 
tinuous tapped or shaken systems, where kinetic energy 
is important. Mechanical spectroscopy using a torsional 
oscillator on granular packings fluidized by continuous 
vibrations have found a growing relaxation time as the 
vibrational amplitude is turned down [72] [73] [74] . It will 
be interesting to compare the rotational and volumetric 
susceptibility, which one can take provisionally to be the 
analogs of viscosity and specific heat, in the same system. 
Do the timescales for the two processes track each other 
as the vibrations are turned down, or do they decouple 
instead? 

What do particle rearrangements in granular com- 
paction look like? In particular, is the ansatz that low- 
frequency responses are contributed by large-scale parti- 
cle rearrangements correct? Visualizing particle motions 
in a 3D packing over the course of an extended experi- 
ment will be very challenging. One possibility is to do 
so on a coarse-grained way using clear particles mixed 
with a small fraction of opaque tracers. By imaging the 
pack using backlighting from multiple angles, a rough op- 
tical "tomography" may be constructed, and collective 
particle rearrangements may be inferred from frequency- 
locked cross-correlations in the image field. Larger sys- 
tem with larger D/d should be studied to see when do 
finite size effects disappear completely. This will give us 
an upper critical length scale for collective rearrangement 
in granular compaction. It will be especially interesting 
if this upper critical length turns out to be very large. 
Since granular compaction is arguably closer to the ideal 
T = jamming transition than either colloidal glasses 
or glass-forming liquids, it is possible that a diverging 



length scale associated with point J will leave a larger 
vestige in granular compaction. 

Throughout this paper I have made purely operational 
analogies between granular compaction and the dynam- 
ics of structural and spin glasses. But there is no first- 
principle justification that the athermal, dissipative pack- 
ing of grains can be described by statistical mechanics 
using macrosocopic observables. That the relationship 
between the volumetric susceptibility and steady-state 
volume fluctuations appear to obey an effective FDT is 
consistent with Edward's proposal for granular statistical 
mechanics [75]. The validity of the statistical mechanics 
approach can be tested in detail if "susceptibilities" in ad- 
dition to Xv can be proposed and measured. Will these 
also obey an effective FDT with respect to the relevant 
fluctuations? If so, is the same effective temperature felt 
for the different susceptibilities? And how does the effec- 
tive temperature depend on the strength and the manner 
in which the packing is driven? It is useful here to recall 
that for all but the simplest thermal systems, the basic 
assumptions of standard statistical mechanics cannot be 
justified except by how well the theory describes experi- 
mental results. Apart from what they can tell us about 
glassy dynamics, one value in experiments such as the one 
here is in providing new and varied data that can test, 
experimentally, the validity of the statistical mechanics 
approach in describing granular materials. 
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